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Abstract. 

We obtain long series expansions for the bulk, surface and corner free energies 
for several two-dimensional statistical models, by combining Enting's finite lattice 
method (ELM) with exact transfer matrix enumerations. The models encompass all 
integrable curves of the Q-state Potts model on the square and triangular lattices, 
including the antiferromagnetic transition curves and the Ising model (Q = 2) 
at temperature T, as well as a fully-packed 0(n) type loop model on the square 
lattice. The expansions are around the trivial fixed points at infinite Q, n or 1/T. 

By using a carefully chosen expansion parameter, g <C 1, all expansions turn 
out to be of the form nfcLi(-'^ ^ gfc)"fc+fc/3fc^ where the coefficients ckfc and (3k are 
periodic functions of k. Thanks to this periodicity property we can conjecture the 
form of the expansions to all orders (except in a few cases where the periodicity is 
too large). These expressions are then valid for all < q < 1. 

We analyse in detail the g — > 1^ limit in which the models become critical. In 
this limit the divergence of the corner free energy defines a universal term which 
can be compared with the conformal field theory (CFT) predictions of Cardy and 
Peschel. This allows us to deduce the asymptotic expressions for the correlation 
length in several cases. 

Finally we work out the ELM formulae for the case where some of the system's 
boundaries are endowed with particular (non-free) boundary conditions. We apply 
this in particular to the square-lattice Potts model with Jacobsen-Saleur boundary 
conditions, conjecturing the expansions of the surface and corner free energies to 
arbitrary order for any integer value of the boundary interaction parameter r. 
These results are in turn compared with CFT predictions. 
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1. Introduction 

Over the years much effort has been devoted to the study of boundary effects in 
statistical models. In the particular case of critical two-dimensional systems a huge 
amount of knowledge has been obtained by the application of the powerful techniques 
of integrability and conformal field theory (CFT) [Il[2]. The critical fluctuations near 
a boundary have been shown to define various boundary critical exponents, many of 
which can be computed exactly. 

Once a model has been shown to be exactly solvable (integrable), it is usually 
rather simple to obtain the dominant asymptotics for the partition function, i.e., 
the bulk free energy /b. The usual route is to diagonalise the transfer matrix in a 
cylinder geometry, using the Bethe Ansatz, and to extract its dominant eigenvalue in 
the thermodynamic limit. Similarly, if the boundary conditions are compatible with 
integrability (in the sense of Sklyanin's reflection equation |3]) one can diagonalise 
the transfer matrix in a strip geometry and retrieve the surface free energy /g. 
Obviously /b and fs give important information about the particular lattice model 
being studied. However, both quantities are non-universal and as such teach us 
nothing about the critical system that emerges in the continuum limit. 

If the lattice model under study is defined on a large rectangle, the next 
subdominant contribution to the free energy comes from the corners. Cardy and 
Peschel [Ij have shown that this corner free energy fc is of universal nature. Indeed 
it contains information about the central charge c of the field theory that governs 
the continuum limit. 

Unfortunately only very little is known about fc from the point of view of exactly 
solvable models. To compute it from the Bethe Ansatz would require knowledge 
about how to implement a given boundary condition in terms of Bethe states. 
Typically this would call for information about all the eigenstates, not just the one 
defining the leading eigenvalue. This obviously defines a very hard problem which, 
as far as we know, has not yet been addressed. 

One of the authors has recently conjectured exact product formulae for e-^*", 
e-^" and e-^'^ in a particular two-dimensional model, viz. the zero-temperature 
antiferromagnetic Q-state Potts model on the triangular lattice [5]. To obtain 
these expressions, the first step was to compute the first 40 terms of their large-Q 
expansion, by combining Enting's finite lattice method (FLM) [6] with exact transfer 
matrix enumerations. The next step was to notice that when reexpressed in terms 
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of the variable x <^ 1 defined hyQ = 2 — x — x ^ these series could be rewritten as 
product formulae of the type 

oo 

llii-x'^r, (1) 

k=l 

where crucially the exponents turned out to be periodic functions of k. Since the 
observed periods (6 for /b, and 12 for /g and fc) were much shorter than the number 
of available terms, it then became feasible to conjecture the exact product formulae, 
valid for any x in the range < x < 1. The model is non-critical in that range, but 
goes to a critical theory when x — 1^ (i.e., Q — i- 4+). 

The purpose of this paper is to extend this type of results to a whole range 
of two-dimensional lattice models. We recall the definitions of the models to be 
studied in section [2] below. They encompass all integrable curves of the Q-state 
Potts model on the square and triangular lattices, including the antiferromagnetic 
transition curves and the Ising model {Q = 2) at temperature T, as well as a fully- 
packed 0(n) type loop model on the square lattice. The expansions are around the 
trivial fixed points at infinite Q, n or 1/T. 

A crucial ingredient in carrying out this programme is obviously to identify the 
correct expansion variable x (that we shall call q in the general case). We wish q to 
have the property that the model is trivial for q = 0, non-critical for < g < 1, and 
critical in the limit g — ?■ 1~. Below we shall focus only on models which are known 
to be integrable (at least in the bulk), and we can therefore take inspiration for the 
choice of q from the exact solution (for the bulk properties). Incidentally, we have no 
reason to believe that nice product formulae hold for models that are not integrable. 

In all cases that we have investigated, it appears that e-^'^, e^^'' and e-^'^ can be cast 
as exact product formulae of the more general form 

oo 

llil-q'T', with ak = Pkk + Ik, (2) 

k=l 

where both sets of coefficients Pk and are periodic functions of k with the same 
period. We find for all Potts and fully-packed loop models included in this study 
that the period associated with e-^" and e-^'^ is precisely twice the period associated 
with e'^'". For Ising models we find on the contrary that all three quantities have the 
same period. At present we have no satisfactory explanation for this observation of 
period doubling for the boundary related properties. 
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Just like in [5] we can then conjecture the exact product formula, provided that 
the series obtained from the FLM has sufficiently many terms to cover at least one 
period (plus a few extra terms to verify the assumption of periodicity, and to account 
for some simple extra factors in front of the product which are sometimes present). 
Only in a few exceptional cases do the series turn out to be too short, but we can 
then at least state how many further terms would be needed in order to establish 
the product form. 

The organisation of the paper is as follows. In section [2] we precisely define the 
models to be studied in this paper. In section [3] we first review the FLM for models 
with free boundary conditions. Then we present our main results in the form of 
product formulae for the bulk, surface and corner free energies of the various models. 
We give details about the period doubling phenomenon for the boundary related 
properties. When studying the critical limit g — !■ 1^ we pay special attention to the 
divergence of the corner free energy, since this can be compared with CFT results [1] . 
As a by-product we obtain results about the asymptotic behaviour of the correlation 
length. 

In section |4] we study the effect of imposing particular (non-free) boundary 
conditions on some of the system's boundaries. We adapt the FLM formalism to this 
case and separate the contributions from corners of different types, i.e., where two free 
(resp. two particular, resp. one free and one particular) boundary conditions meet. It 
is shown analytically that the contributions to the free energy from corners of different 
types simply add up, as expected for an inherently non-critical system. We apply 
this formalism to a family of Potts-model boundary conditions recently introduced by 
Jacobsen and Saleur (JS) [7] in which a parameter r controls the weight of Fortuin- 
Kasteleyn clusters that touch the particularised boundaries. Explicit results are 
obtained for any integer value of r. To finish, we comment on the relation with CFT 
results for such boundary conditions. 

2. Models 

In this section we define the models to be studied in this paper and briefly review 
their most relevant properties. These models constitute all integrable cases of the 
Potts model on the square and triangular lattices. We also study a model of two 
mutually excluding sets of fully-packed loops, known as the FPL^ model |l8j. 

In the entire paper we will define the (dimensionless) free energy as the logarithm 
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of the partition function, / = InZ, i.e., without the conventional minus sign. 
2.1. Potts model: generalities 

Let us recall that the Q-state Potts model is a model of interacting spins on a lattice, 
allowing each spin to be in one among Q different states, and such that the interaction 
between neighbouring spins depends only on whether they are in the same or different 
states. The underlying lattice, or graph, is denoted as G = {V,E), where V (resp. 
E) is the set of vertices (resp. edges). The associated (dimensionless) hamiltonian 
thus reads 

(3H = -K J2 a, = l,...,Q, (3) 

where K = J/ksT is the (dimensionless) coupling constant. 

It is well-known that the partition function can be rewritten as [9] 

z=J2 Q'^^^^'^' , (4) 

ACE 

where the sum runs over subsets A (of cardinality \A\) of the lattice edges E. Each 
connected component in A is known as a Fortuin-Kasteleyn (FK) cluster, and k{A) 
denotes the number of connected components. The temperature parameter v is 
defined as f = e-'^ — 1. 

Obviously, the FK formulation makes sense also when Q is an arbitrary real 
number. Suppose now that the Potts model is defined on an infinite regular lattice 
(|l^|, \E\ — 7- oo) with coordination number z = 2\E\/\V\. If we study the model on 
some curve in the {Q,v) plane with asymptotics 

v^Q^/' for|g|>l, (5) 

then there will be precisely two dominant contributions to Z: A = E (with k{A) = 1) 
and A = (with k{A) = \V\). In other words, there is phase coexistence between the 
low-temperature and high-temperature phases. All subdominant contributions can 
be written perturbatively as powers of the small parameter 1/Q. This perturbative 
picture is the starting point for making series expansions of the free energy. 
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Figure 1. Different types of lattices to be used througliout this work. From left 
to right: a) square lattice in a rectangle of size M x N, h) triangular lattice in 
a (deformed) rectangle of size M x N, and c) triangular lattice in an equilateral 
triangle of size M. 

2.2. Lattices and their orientation 

In a detailed study of boundary effects it is important to specify not only on which 
regular lattice the model is defined, but also how the boundaries are oriented with 
respect to the lattice's symmetry axes. 

Moreover, we do not expect that the free energy series can be written as nice 
product formulae of the type ([2]) for an arbitrary lattice model of interest. The 
model will in general have to be integrable to yield nice expressions. We recall 
that a model remains integrable in the presence of boundaries, only if the boundary 
conditions satisfy the reflection equation [3]. This places strong constraints on how 
the boundaries can be oriented. 

The Potts model on the square and triangular lattices is integrable for certain 
curves in the {Q, v) plane, which we briefly review below. Boundary integrability 
turns out to hold as well, provided that the boundaries are oriented parallel to the 
principal axes of the lattice. In particular, an M x rectangular piece of the square 
lattice will be oriented as shown in Fig. [1^. 

When dealing with a triangular lattice, we can again take an M x (deformed) 
rectangular piece, as shown in Fig. [lb. This lattice can be considered simply as a 
square lattice with added diagonals. This point of view is often convenient, since then 
the FLM formulae can be taken over from the square-lattice case, and the transfer 
matrix algorithms need only very minor modifications. 

There is however one disadvantage: since the four corners are not equivalent, the 
corner free energy will be a mixture of two contributions from corners sustaining an 
angle 7r/3 and two contributions from 2tt/3 corners. To separate these contributions 




Figure 2. Phase diagrams in the {Q, v) plane of the Potts model on a) the square 
lattice, and b) the triangular lattice. The models are integrable along the curves 
shown. Critical (resp. non-critical) behaviour is signalled by solid (resp. dashed) 
line style. For further details on these curves, their properties and nomenclature, 
and the renormalisation group flows, the reader can consult [10l[Il] and references 
therein. 



we shall find it advantageous to consider as well the triangular lattice inscribed in 
an equilateral triangle of side length M, as shown in Fig. This involves three 
7r/3 corners, but makes both the FLM formulae and the transfer matrix algorithms 
slightly less performing. 

2.3. Square-lattice Potts model 

The phase diagram of the square-lattice Potts model is shown in Fig. [2^. The model 
is integrable [I2] along the selfdual curves 

v = ±^/Q, (6) 

corresponding to a critical (resp. non-critical) theory for < Q < A (resp. Q > A). 
The conformal properties in the critical regime are known [13]. The only simple 
property that we need here is that the point {Q, v) = (4, 2) is critical with central 
charge c = 1. 

Note that satisfies the criterion ([5]), thus making possible a perturbative 
expansion for Q ^ 1. We shall see below that this expansion enables us to attain 
the point {Q,v) = (4,2) as an appropriate limit. 
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The square-lattice Potts model is also integrable [H] along the mutually dual 
antiferromagnetic curves 

v = -2±^A-Q, (7) 

corresponding to a critical (resp. non-critical) theory for < Q < 4 (resp. Q < 0). 

The conformal properties along ([7]) are quite intricate, but much is known 
[ini [m [151 US]. In particular one has c = -1 in the limit (Q, v) (0, 0). Since (JZD 
again satisfies ([5]) the perturbative expansion enables us to attain this point as an 
appropriate limit. 

Note finally that the Q-state Potts model is equivalent to a problem of self- 
avoiding loops [17] . The loops are defined on the medial lattice and each one carries 
a weight n = y/Q. The loop formulation can further be brought in equivalence with 
the six- vertex model [TTj. The six- vertex model is homogeneous on the curve and 
staggered on ([7]) [18]. 

2.4- Triangular-lattice Potts model 

The phase diagram of the triangular-lattice Potts model is shown in Fig. The 
model is integrable along the cubic curve [19] 

v^ + 3v^ = Q. (8) 

Once again is satisfied. 

For < Q < 4 the critical behaviour for v > (i.e., along the upper branch of 
(E])) coincides with that of the square lattice. Along the lower branch of ([H]) one finds 
[TT] the same critical behaviour as on the antiferromagnetic curve ([7]) on the square 
lattice. 

The two critical points that can be accessed perturbatively are {Q,v) = (4, 1) 
with c = 1, and {Q,v) (0, -2) with c = -1 [IT]. 

We shall sometimes refer to (|H1) as the selfdual curve. Indeed, if one combines a 
duality transformation of the triangular-lattice Potts model with a partial summation 
(decimation) over one half of the spins in the dual model on the hexagonal lattice, 
the result is a non-trivial transformation whose fixed points are exactly ([8]). In this 
sense the nomenclature "selfdual" is justified when speaking about the curve ([8]). 

In the special case v = —1 and any real Q, the Potts partition function Z 
reduces to the so-called chromatic polynomial. For integer Q this can be interpreted 
as a colouring problem. To be precise, Z equals the number of proper vertex 
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colourings of the lattice, i.e., assignations to each of the vertices in G of one among 
the Q different colours, in such a way that adjacent vertices carry different colours. 

The chromatic polynomial is integrable [201 121] for any real Q. It corresponds 
to a critical theory if and only if < Q < 4. Its phase diagram within the critical 
phase is complicated [21] and has been reviewed in [22]. The only properties that 
we shall need here are that the point [Q, v) = (4, —1) has c = 2, while the limit 
{Q,v) —J- (0, —1) corresponds to c = —1. 

The bulk, surface and corner free energies for the chromatic polynomial on the 
triangular lattice have already been discussed in [5]. In section l3. 2. 4l we shall however 
be able to go beyond these results and separate the contributions to the corner free 
energy for two different types of corners (of angles ^ and | ) . 

2.5. Ising model 

For Q = 2 the Potts model reduces to the Ising model with (dimensionless) 
hamiltonian 

I3H = -Ki,^^ S,S,, S, = ±l. (9) 

The coupling constants in ([3]) and @ are obviously related by -ft'ismg = ^K. 

The Ising model is integrable at any temperature on both the square and the 
triangular lattices [18]. We shall only be interested in the cases -Raising > 0. The 
critical point on the ferromagnetic critical curves is situated at -ft'ising = |log(l + A/2) 
on the square lattice (see Fig. |2^) and at -ft'ismg = \ log 3 on the triangular lattice 
(see Fig. |2]d). The low-temperature perturbative expansion gives access to the critical 
point in an appropriate limit. 

When studying the Ising model, we will always treat it as a Potts model and 
rewrite -falsing in terms of K. 

2.6. FPL^ model 

Apart from the special cases of Potts models described above, we also consider a 
model of two mutually excluding sets of fully-packed loops on the square lattice, 
known as the FPL^ model [8]. The allowed configurations are such that each vertex 
is in one of the six states shown in Fig. [31 A weight ni (resp. 712) is associated with 
each black (resp. red) loop. 
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Figure 3. Allowed vertex states in the FPL^ model. 
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Figure 4. A finite bipartite 4-regular lattice of width M = 3 and height N = 2. 

The FPL^ model is integrable when the two weights are equal, ui = n2 = n 
[23| [211 125]. We study henceforth only this integrable case. It gives rise to a critical 
(resp. non-critical) model when —2 < n <2 (resp. \n\ > 2). 

When the FPL^ model is defined on a bipartite 4-regular lattice, its 
configurations are in bijection with those of a 3-dimensional interface model on the 
dual lattice |8]. This allows for the redistribution of the loop weights ni and n2 as 
local complex vertex weights. In the critical regime the resulting local formulation is 
the key to the exact derivation of the critical exponents by the Coulomb gas method 
[261 ET]. 

To investigate the surface properties of the FPL^ model without breaking any 
symmetries (which might lead to a change of universality class) one therefore needs 
to define it on a finite lattice which is bipartite and 4-regular [28]. An appropriate 
choice of lattice is shown in Fig. HI 

3. Bulk, surface and corner free energies with free boundary conditions 

In this section, we compute the bulk, surface and corner free energies of each of the 
considered models in the case of lattices with free boundary conditions. Achieving 
the goals described in the Introduction relies in all cases on writing these free energies 
as power series in terms of some good perturbative parameter. Long series are then 
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Figure 5. Representation of the ground state used as a basis for the finite lattice 
method (left), and of one of the first excitations to be taken into account in the 
series expansion (right), in the case of the square-lattice Potts model. The black 
edges on the right correspond to the subset ACE appearing in The blue 
box surrounding A corresponds to the smallest finite sublattice in which this 
excitation appears. 



generated by combining transfer matrix techniques with the use of Enting's finite 
lattice method [6], which we first review below. 

For each model being considered, the first term of the series expansion is some 
unique (or sometimes twofold degenerate) ground state. As mentioned at the end of 
section I2T3] we shall treat the Potts models dl]) in terms of the equivalent loop models 
on the medial lattice. The ground state is then chosen as the unique state in which 
one loop surrounds each vertex in the original lattice G; see Fig. [51 Equivalently, the 
edge subset ACE appearing in (jl]) is empty, A = (/). 

Similarly, for the Ising model the ground state is chosen as one of the two 
equivalent ferromagnetic states with all spins pointing in the same direction. 

For the ferromagnetic Potts model (both on the square and triangular lattices), 
the chromatic polynomial on the triangular lattice, or for the FPL^ model, a good 
expansion parameter will turn out to be, rather than the inverse of the loop 
fugacity {n = y/Q), the parameter q defined by 



Q = q + -. 



(10) 



Actually, for these models q is nothing but the quantum group deformation 
parameter. Indeed, all these Potts models present a quantum group symmetry 
Uq{sl2)- The FPL^ model is endowed with the larger symmetry Uqi^sl^). 
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For the critical antiferromagnetic Potts model on the square lattice the suitable 
parameter q is defined by 



Finally, for the Ising model the choice of the expansion parameter q is far less 
obvious. Indeed, the approach outlined in the Introduction makes it necessary that 
the critical point correspond to the limit q ^ 1^. The correct choice of q then 
involves parameterising the coupling constant in terms of suitable elliptic functions, 
as in Baxter's analytical results [291 HE]. We defer further details to the corresponding 
sections below. 

3.1. FLM for lattices with free boundary conditions 

We now outline the principles of the FLM for lattices with free boundary conditions. 
The formalism depends on the choice of regular polygons in which the excitations 
over the ground state are inscribed. The most common and well-known choice 
is rectangles, as illustrated by the blue box in Fig. |5l We review this first (in 
section lS.l.ip for convenience, and in order to fix the notation. Then, in section [3.1.2[ 
we derive the required modifications for the case of triangular lattices with excitations 
inscribed in equilateral triangles, cf. Fig. [T]:. 

3.1.1. Lattices inscribed in a rectangle. We here give a brief review of the earlier 
work by Enting [6] concerning the basis of the finite lattice method. The FLM 
allows one to compute the free energies of statistical models defined on an infinite 
rectangular lattice as a series expansion in terms of finite rectangular graphs. Apart 
from the square lattice, this method can also be applied to the triangular lattice, 
which is then considered as a square lattice with added diagonals; see Fig. [Tb. 

In practice, we will make extensive use of the FLM to approximate the bulk, 
surface and corner free energies of an infinite M x N lattice, such as depicted in 
Fig. [T^-b, from free energies for finite m x n lattices. The latter are computed by 
using numerical (but exact) transfer matrix enumerations. The bigger the m x n 
lattices whose partition function can be computed numerically, the longer will be the 
resulting FLM series. 




(11) 
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The basic idea underlying the FLM is to write the free energy for the M x N 
lattice as a sum over all possible sublattices, 

fM,N= Yl E (M-^ + l)(iV-j + l)/,,, (12) 

[i,j]c[M,N] i<M,j<N 

where the notation [i,j] C [M,N] means that the summation is performed over 
all sublattices of size i x j that lie inside the M x N lattice. The contributions 
fij — that must not be confused with the lattices' free energies fij — are defined 
self-consistently as functions of the latter by inverting the equation 

fm,n = fi,j = ("^ - ^ + - J + • (13) 

[i,j]C[m,n] i<m,j<n 

This results in 

Ai = fm,nV{m,i)r]{nJs) , (14) 

m<i,n<j 

where the functions r] are defined by Enting [B] as 

{1 if m = i or m + 2 = i and i > 2 , 
2 if m + 1 = z and z > 1 , (15) 
otherwise. 

For an infinite M x N lattice, the sum f ll2p is restricted to a certain range of 
finite mxn sublattices by imposing some cutoff set. In the course of this paper, and 
as prescribed in [6], this cutoff set is chosen to be 

B{k) = {[m,n] ,m + n = k} . (16) 

Hence, (fT2|) formally reads 

fM,N^ Y {M-z + l){N-j + l)f\,, (17) 

[i,j]<B{k) 

which is generally shown to give an approximation of fM,N up to an order increasing 
with k in terms of some good perturbative parameter, as discussed above and detailed 
in the following sections. 

The free energy (flTl) can be written as the sum of a bulk contribution MNf^, a 
surface contribution (M + N)fs, and a corner contribution /c, that is 

fM,N = MNf^ + (M + N)f, + /, (18) 
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These three contributions can be computed separately. The final form of the FLM 
expansions with cutoff B{k) come out as [6j 

fh ^ ^ fm,n(^^m,k—n 3(5.m,fc— n— 1 ~l~ '^^m,k—n—2 ^m,fc— n— s) ; 

[m,n]<B(k) 

fs = /m,n((l — ^)^m,k-n + (3m — l)6m,k-n-l , 

[m,n]<B{k) 

- (3m + l)6m,k-n^2 + (m + l)6m,k~n-3) (19) 

fc= ^ /m,„((m - l)(n - l)(5m,fc-n + (1 + m + n - 3m?2)(5m,fc_„_i 

[m,ra]<B(A;) 

+ {3mn + m + n- l)Sm,k^n-2 - (m + l)(n + l)5m,k-n-3) , 

where all the fm,n with m or n non-positive are zero by convention. Compared to 
Enting's formulae, the expression for fs has been simplified under the assumption 
that the fm,n are symmetric in m and n, which will always be the case in this paper. 

3.1.2. Triangular lattices in triangles. For the sake of studying models defined on 
the triangular lattice, it is sometimes useful to work with only one type of corners. 
One then seeks the free energy fM of the model defined on a large equilateral triangle 
of side M, as shown in Fig. [It. We have derived the corresponding FLM formulae, 
going through the same steps as above. 
The starting point is now 

i<k 

where the self-consistency condition 
yields 

fi = fm{^m,i ~ 35m,i-l + 35m,i-2 ~ ^m.j-s) • (22) 

m<i 

Finally, the free energy is written as 

/M = ^^^^^V b + M/3 + /c, (23) 
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where the different contributions are given by 




m<k 



/s = ^ /m((l - rn)5m,k + (2m + l)5m,k-i - (m + 2)5m,k-2) 



(24) 



m<k 



/c = 2 X] fm{{m'^ - 3m + 2)6^^k + 2(1 - m^)6m,k-i 



m<k 



+ (m^ + 3m + 2)5^^fc_2) . 



3.2. Results 

3.2.1. Selfdual Potts model on the square lattice. On its selfdual curve, the square- 
lattice Potts model is equivalent to a loop model, where each loop carries a statistical 
weight n = \fQ. The partition functions of finite lattices of size M x N reads 



where Pd{x) is a polynomial of degree d = MN—1. The constant coefficient -Prf(O) = 1 
corresponds to the ground state shown in Fig. [5^. 

The FLM formulae thus reduce to an expansion in powers of 1/n, which is 
correct up to an order determined by the size of the graph with highest weight that 
cannot fit into the cutoff set B{k). Compared to the ground state configuration, this 
property is satisfied by the graphs in which the edge subset A C £" of (jlj) forms one 
straight or L-shaped figure of total length k + 1, which contribute to the order n"*^. 
We thus expect the FLM formulae f|T9|) with a cutoff B{k) to give an approximation 
of the M X N free energy that is correct up to order n~^. 

Using a transfer matrix algorithm, we have computed the polynomials Pd{x) 
for all finite lattices within the cutoff set 5(31), allowing correct results up to order 

Using instead of 1 /n the small parameter q defined by 



the exponentiated free energies e-^^", e-^" and e^^'^ finally turn out to have up to this 
order the following nice expressions 



Zm,n = n^'^'Pdin-') , 



(25) 



1 

n = q-\ — 



(26) 
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1 



1 



fc=i 



8fc-l 



q 



8k-5 



16 



(27) 



k=l 



^8fc-6)(l_^8fc-4)4(l_g8fc-2) 



The expression for e-^"" is shown in Appendix A to be equivalent to the analytical 
expression given by Baxter in [TS], that is 



-kX 



k=l 



sinh 2k P 
cosh kX 



(28) 



where ip = —^^^n-^ooJi^^Zn is defined as the dimensionless free energy in the 
thermodynamic limit (here Zj^r denotes the partition function of a system with 
spins). The correspondence with our notation is 



- fh , 

2 cosh A , 
sinh P 



gV2 sinh(/3 - A) ■ 

The regular form exhibited by the product expression for e-^"" — which is proved 
starting from Baxter's result (128|) in Appendix A.l — and its structural similarity with 
the corresponding conjectured expressions for e-^= and e-'^^ encourages us to assume 
that the two latter expressions are also valid to all orders. These expressions are, as 
far as we know, new|il 

In general we expect e-^*", e-^" and e-^"" for all the models studied in this paper 
to have an expression of the form ([2]), modulo some simple overall factors such as 
gi\q^i)2 {1 — q) appearing in fl27|) . Moreover, we expect the coefficients Pk and 7fc 
to be periodic functions of k. For instance the results found in fl27j) correspond to 
the simplest case of Pk = 0; and the periodicities are 4 for e-^'^, and 8 for e-^'^ and e-^'^. 
We shall discuss the issue of periodicity further in section 13.2.81 

The above assumption will be further corroborated by the study of other models 
in the remainder of the paper. 

We remark that a result such as e-^'' in fl27|) can be said to be established beyond 
any reasonable doubt, since the coefficients 70, 71, . . . , 731 determined numerically 
cover four complete periods. Once the periodicity property is established (or 

I But in the case of e^= it seems likely that an equivalent result exists somewhere in the literature. 
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assumed) the determination of the coefficients over a single period (or a httle more 
to be sure) would be sufficient. 

One more important feature of the generic form ([2]) is that it gives a particular 
meaning to the limit g — )■ 1~, which turns out to be critical for all the considered 
models. The convergence of products of the type ([2]) when q 1~ indeed depends in 
a non-obvious way on the coefficients ak, giving rise to different critical behaviours 
to be studied in the following. 



3.2.2. Antiferromagnetic Potts model on the square lattice. The antiferromagnetic 
transition curve [H] for the square-lattice Potts model is given by ([7]). To resolve 
the square root we parameterise Q = —{q^^ — q)^, with q taking real (resp. complex) 
values in the non-critical (resp. critical) regime. We have then v = (g^^^ — g~^/^)^. 
This allows us to compute the partition function of finite lattices as polynomials in 
q. More precisely, the partition functions for lattices of size M x N read 

M+N-3MN 

Zmn = ^zr^^P^i^) ' (29) 

where Pd{(l) is a polynomial of degree d = 3MN—M—N in q with constant coefficient 
-fd(O) = 1. To obtain P^, FLM calculations are applied in exactly the same way as 
for the selfdual curve. The numerical calculations are now done up to order g^^, i.e., 
using the finite partition functions with the cutoff set -8(23). This gives for the bulk, 
surface and corner and free energies: 



g/b 



g/s 



g/c 



^ k=l ^ ^ 



°° /I „16fc-14Wi „16fc-2\ °° /I „8fc-3\2 

n (i-g ){l-q ) yr / i q \ 



(1 - ^) 

^ (l_g4fc-2)(^_g8fc-4)3 



The presence of an overall minus sign in e-^*" stems from the fact that Q < 0, 
and should not be taken into account as for the critical properties of the model when 
Q ^ 0. The expression found for e-^*" is shown in Appendix A. 2 to be equivalent to 
an expression given by Baxter in section 12.5 of [H], up to an analytic continuation. 
More precisely, using the same notations as for the self-dual case, the free energy ip 
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is given for Q > 4 by fl28l) . where this time (3 has to be replaced by a parameter u 
defined by 

sinhfA + u) ^ , , 

^ =—rr, 0<3^i<7r, 31 

smh(A — u) 

the expression being vahd under the condition < < A. On the contrary, the 
expressions for e-^= and e-^'^ are new. 

3.2.3. Self dual Potts model on the triangular lattice. The selfdual transition curve 
on the triangular lattice is given by the cubic (|8]). It can be resolved through the 
parameterisation 

V = - 1 + t + r^ (32) 

The FLM calculations are then conducted in terms of the variable t = cj^l"^ . 
The partition functions read 

where Pd{t) is a polynomial of degree d = 8{MN — 1) in t with constant coefficient 
-Pd(O) = 1- Finite lattice calculations similar to those used so far are then applied 
to obtain the Pa- The only difference is that this time the highest order graph that 
cannot fit into the cutoff set B{k) is the straight line graph oriented diagonally, 
formed by ~ k/2 diagonal edges, which contributes to order t'' in the partition 
function. We write the final results in terms of the same variable q as used in the 
square-lattice case. 

The numerical calculations are conducted up to order t"^^ (i.e., with cutoff set 
5(23)), that is to order q^^^^. They lead to the following product expressions for the 
bulk and surface free energies: 

g/b = 




= I I I r 1 • (34) 
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The expression (!34|) of e-^*" is shown in Appendix A. 3 to be equivalent to an 
expression given by Baxter in section 12.6 of [ISj, that is, 



_;j;^sinh 2k /3 
k cosh kX 

k=l 



(35) 



with the same notations as in section I3.2.1I The expression found for e-^" seems, on 
the contrary, new. 

The series expansion of e-^'^ appears to be more comphcated. Up to order g^^/^ 
we find: 

= (1 _ _ q^y^l - gi)(l - qf)-\l - q^)-\l - qf)-' 

(1 - gf )-5(l - q^il - gf )-9(l - gf f (1 - q^)-\l - qff 

(1 _ gf - gl0)15(l _ gf _ gf )13(1 _ gl2)~22 

(1 _ )19(1 _ gf _ gl4)19(i _ )-29 ^ ^(^f ) _ (gg^ 

We beheve that the correct expression is still of the form (|2]), but for the first time 

We recall that ( l36l) applies to the triangular lattice inscribed in a rectangle, cf. 
Fig. [T]d. This has two different types of corners: 1) two of angle ^ each connected 
to three edges, and 2) two of angle | each connected to two edges. The contribution 
from each type of corner to e-^'^ is expected to be different. Therefore we should be 
looking for an expansion of the form 

oo 



e 

k=l 



ll{l-qr^'^<\ (37) 



where a^*^ = Pk^k + '-fl^\ In accordance with what was observed for previous models, 
we expect the exponents /3^^ and 7^,*^ to exhibit the same periodicity as in the 
expression for e^^"^; see flMl) . This periodicity is 8 in the g variable, but since powers 
of g^/^ occur we need 12 consecutive terms to determine one period. Moreover there 
are now 4 different types of exponents, so the complete determination would require 
a cutoff set at least -8(48), which is clearly out of reach of the methods used this far. 

To proceed we may turn to lattices inscribed in a triangle [cf. Fig. [1]:], for which 
we expect the bulk and surface free energy to be unchanged, but which possess three 
corners of the same type, namely of angle | (type i = 2). Computing finite-lattice 
partition functions by a transfer matrix algorithm similar to those used previously, 
and using the FLM formulae of section 13.1.21 yields indeed the same expressions as 
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above for e-^'^ and e-^^ up to order t^^, except from the fact that the exponents = 2 
in are now replaced with ak = 3, which stems from the different definitions (JT8l) 
and fl23l) chosen for fs in the rectangle and in the triangle case. For the corner free 
energy we find up to this order 

oo ^ 



k=l 



(2) 

Quite remarkably, this expression shows that = 0, so the linear growth of the 
exponents visible in (!36l) must be due to pj^^ ^ 0. Insight can thus be gained by 
factorising out the contribution of two corners of angle | [that is, eq. f l38p to the 
power |] from (|36|) . The remainder corresponds to the contribution of two corners 
of angle and is given by 

2/3 

-2^ 



4 , 



x(l-gf)-i(l-g2)-i(l-g3) 

.k=l 

(1 - qf^l - - qfni - qf y'^l - q'^^l - qf )-'' 

(1 _ _ ^12)-22(^ _ )21(^ _ _ ^14^19(1 _ qfy29 

+ 0{qf). (39) 

Unfortunately this is not quite enough in order to be able to conjecture the 
corresponding exact expression. According to the above discussion, we would need 
the first 24 exponents (and ideally a few more as a verification), but we only have 22. 
For the sake of checking consistency, we tried to adapt the FLM formulas to lattices 
inscribed in hexagons, that is, lattices with six corners of type i = 1. However, it 
seems that in this case the FLM formulae are impossible to invert. 

3.2.4- Chromatic polynomial on the triangular lattice. The finite partition functions 
for the Q-colour chromatic polynomial on a triangular lattice inscribed in a rectangle 
were calculated by one of the authors in |5]. This resulted in the following product 
formulae for the bulk, surface and corner free energies: 



g/b 



1-Q (l_a;6fc-3)^i_3,6fc-2)2^i_^6fc-l 
k=l 



-in 1 

oo 



e 

k=l 



a;m-9)(i _ a;i2fc-8)(i _ a;i2fc-7)(i _ a;m-2) 
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k=l 

(1 _ a;i2fc-3)5(i _ a;m-2)2^^ _ 3,m-i)2(^ _ ^m^ ^ ^^^^ 

where the expansion parameter x is defined this time by(5 = 2 — x — ^. 

As for the selfdual Potts model on the triangular lattice, these expressions can 
be checked by FLM calculations on the triangular lattice inscribed in a triangle. We 
have computed the corresponding finite-lattice partition functions up to triangles of 
size M = 21, which gives access to the bulk, surface and corner free energy series 
up to order ~ The expressions for e-^'' and e-^= coincide with fHOl) as expected, 
whereas we find for the corner free energy 



e^^ = 11(1 - xi2^-i°)(l - x'^'^-' fil - a;i2'=-«)(l - x'^'^-'f 
k=i 



(1 - x'^'-^)il - x'^'^-' fil - x'^') . (41) 

We can now combine eq. fj40|) for two corners of type i = 1 and two corners of 
type i = 2 with eq. fHTl) for three corners of type i = 2, to obtain the exact formulae 
for a single corner of each type: 



k=l 



(1 _ 3.m-5)-l/2(^ _ ^m-4)-5/6(^ _ ^m-3)3/2 

(l_a;i2/c-2)2/3(^_^my/6^ (42) 

oo 

e/^(f) = J](l-a;i2^-i0)i/3(i_3,m-9)2/3(^_^m-8)i/3 

k=l 

(1 - x'''-'){i - x'^'-^r'/'ii - x^'"-") 

(1 _ 3.m-4)-2/3(^ _ 3,m-3)(^ _ ^m-2)i/3 
(l_a;i2fc-i)(i_3.m)i/3^ (43) 

The chromatic polynomial on the triangular lattice is critical in the region 
Q G [0,4]. As for the selfdual Potts model, we are interested in the limit Q — )■ 4"*", 



Corner free energies and boundary effects 



22 



which we parameterise in terms of the usual variable q. It is related to the above 
variable x hj x = —q^. To turn the above products into expressions in terms of q 
thus amounts simply to the following replacements, 

1 — x^^ is replaced by 1 — , 

1 — x^^^^ is replaced by . (44) 

We will also be interested in the limit Q — ?■ 0^, which can be recovered as g — i x 1~, 
i.e., using complex values of q. 

3.2.5. FPL^ model. For lattices of the type shown in Fig. Hlthe partition function 
is 

Zmn = 2n2^^+*^+^P,(n-2) , (45) 

where Pd{x) is a polynomial of degree d = MN + [(M + N)/2\ — 1 in the variable 
X = with constant coefficient -Pd(O) = 1. The overall factor of 2 in fH5l) is 
unimportant for what follows, and it is convenient to regard the power of n as a 
multiplicative contribution of to e-^'' and a contribution of n to e■^^ 

We have computed the polynomials Pd in (l45ll for all lattices with M + N < 17. 
Introducing the usual parameterisation n = q + q~^, and using the FLM, this yields 
series expansions in powers of q for the the bulk, surface and corner free energies, 
e-^*^, e-^" and e-^'^, that are correct up to order g^^. Since Pd depends on n~^, only even 
powers of q appear in these series expansions. 

Exact product formulae for these series can readily be conjectured. Reinstating 
the above-mentioned multiplicative contributions we obtain 



g/b 



g/s 



9 fciv (1 -g^^-6)(l -g^'=) ; 



^ k=i 



J< 



k=l 



TT / (l-gl6.^1^)(l-gl6.^10)(l-gl6.-2)3 X ^ 
~ n (I _ ql6k-Uy(^l _ gl6k-12-j^l _ gl6fc-10-) ^ 



'I — qiek-Gy^^l _ gl6fc-4^5^]^ _ gl6fc-2^ • 
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Note that the above expression for e-^*^ is identical to the one obtained in Eq. (64) 
of [23] from an exact Bethe Ansatz solution. On the contrary, the expressions found 
for e-^= and e-^'^ are new. 



3.2.6. Ising model on the square lattice. Turning to the case of the two-dimensional 
Ising model on the square lattice with nearest-neighbour coupling J, the partition 
function is calculated from the FLM formulae as a low-temperature expansion in 
powers of the variable x = = e~^'^ . The partition functions read 

^ Pd{x), (47) 



g.MN-M-N 



where Pd{x) is a polynomial of degree d = 2{2MN — M — N) in x with constant 
coefficient -Pd(O) = 1. From the same arguments as for the selfdual Potts model, the 
use of a cutoff set B{k) leads to an approximation that is correct to order ~ x^'^. 
The choice of a "good" variable in terms of which the bulk, surface and corner free 
energies take a well-behaved factorised form is inspired by Baxter, Sykes and Watts 
We therefore set 



oo 



;l-g8fc-7)(l_g8fc-l^ 



-2 ^ ^1/2 ■ 

k=l 

One can verify that the critical point /ScJ = \ log(l + \/2) corresponds to the limit 
g — 7- 1^. In terms of the variable q we find, using the cutoff set 5(49) (i.e., up to 
order g^^), that 



g/b 



g/s 



gl/2 11 ^X — g8fc-7^8fc-7('X _ g8fc-3')8fc-3 — q^^-^){\ — 



g/c 




(49) 

8fc-4 



fe=l 

The expression found for e-^'' is shown in Appendix A. 4 to be equivalent to the 
analytical expression given in section 11.8 of [18] for the dimensionless free energy in 
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the thermodynamical limit, namely 

^ = _2K-T + g{z)+g{z''), (50) 
where K = (5J is the dimensionless Ising coupling, 2; = 1 at the critical point, 

m=l ^ J / \ -1 / 

and r can be expressed function of q as 

m(l — q' , 

m=l ^ ^ ' 

Another equivalent expression was given in [29], namely 



00 



r = y Y~\ . (52) 



n=l 



^-^^+2^ n(l + g2n)(l_^8n^ " ^^"^^ 



3.2.7. Ising model on the triangular lattice. In the case of the Ising model on a 
triangular lattice, the natural variable to be used for product expansions is inspired 
by section 11.8 of Baxter's book [18], namely 



.^^,m(^-^"-":)(^-^""-^ (54) 

The critical point Kc = (3cJ = log 3 corresponds once again to the limit g — t- 1^. 

The finite-lattice partition functions read 

2 

^^^^ = 3,3Af7V-2Af-3iV+l^'^(^^) ' (^^) 

where Pd(x^) is a polynomial of degree d = 2MN — M — N m with constant 
coefficient -Pd(O) = 1. From the same arguments as for the triangular- lattice Potts 
model, the highest order graph that cannot fit into the FLM cutoff set B{k) is 
a diagonally oriented line graph of length ~ A;/2, corresponding to inverting ~ k 
surrounding links and thus contributing to the partition function up to an order 
~ x^'^. Leading calculations up to the cutoff set -8(53), we thus find up to order g^o/^ 
for the bulk free energy 

6fc-l 



k=l 



(1 -g^'=-3)(l-g«'=+3)(l-g 




n ^n^^ I n 

11(1_^8.-6)(1_^8.-2)11 
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This expression is shown in Appendix A. 5 to be equivalent to the analytical 
expression given in [18] for the dimensionless free energy in the thermodynamic 
limit, namely 

ip = -3K - T + 3g{z) , (57) 

with the same notations as on the square lattice, and this time z = q~^/^ at the 
critical point. Another equivalent analytical expression was given by Baxter, Sykes 
and Watts in [29]. 

We were unable to obtain a regular factorised form for the surface and corner 
energies, for the same reasons as those given in section 13.2.31 in the case of the 
triangular-lattice Potts model. We nevertheless give hereafter the first terms of the 
factorised series expansion in terms of the variable q: 



e^= = x\l - q^y\l - qi)\l - - g^)-^(l - g^^"^ 



3 

;i _ )2(1 _ qff^l _ gf )-6(l _ gf )-6(l _ )12 
;i _ gf )8(1 _ gf )-10(l _ gf )-12(l _ )6(i _ )12 ^ 



ef^ = -(1 - - _ _ qh-\l - q 

X 

''^-4/1 ^l\4/i '3\-2n ^^\3n ^4n-3 



1 - q-3)-\l - q^ni - q-T - - q' 

1 _ qf)-\l - _ qf)-\l - gir)-^(l - q')^ 

— ^-4/1 ^W-\lOn ^7\-2/i ^^\6/i ^8\-2 



3 



(1 - - qr'a - - - 
(1 - q'^r\i - q'-^r'ii - q'r'a - q^Ya - - ^ 
(1 - Q'r'a - - - Q'ra - ^^^r 

As for the previous models defined on triangular lattices, FLM calculations 
were also performed on lattices inscribed in triangles. Having computed finite-lattice 
partition functions for triangles of size M < 23, which yield series expansion up to 
order x^^, that is q^ , we find for e^^'^ the same expression as fl5^ . The first terms 
of e-^= agree with (158|) . except for the fact that all exponents are here multiplied by 
I due to the different definitions (fT8|) and (123!) chosen for /g. Concerning the corner 
free energy, we find in this case 

e^- = 2(1 - g^)~^(l - gi)^(l - q^f{l - q^y\l - q^y^ 

(l-g3) _ g 3 )y(l - g 3 )y(l - g 3 ) ^0(l_g3) 
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(1 _ )21(1 _ )15(l _ qf)-\l _ gf _ gf )-24 

+ 0{qf). (59) 

To fix the corresponding factorised form, we first note that eqs. fj49l) and fl56|) giving 
e'^'" for the Ising model on respectively the square and triangular lattices are both 
expressions of the general form ([2]) with exponents and 7^ that are 8-periodic in 
q. Seeing that e-^'^ for the square-lattice Ising model is also of this form, still with 
periodicity 8, we can assume the same to hold on the triangular lattice. The above 
development fl59l) falls slightly short of furnishing enough coefficients in order to fix 
Pk and 7fc unambiguously. However, it is consistent with the following appealing 
conjecture 



12fc-9 / „, 9 \ 12fc-3 



g/c 




(60) 



Our confidence that ( I6OI) is indeed correct in enhanced by the fact that its asymptotic 
behaviour predicts the same universal divergence of the correlation length on the 
square and triangular lattices (see section [3.4.5p . 

3.2.8. Periodicities of product exponents. Table [T] reviews the periodicities in terms 
of the variable q that we have observed in the exponents for the various models 
described above. It is seen that for the Potts and loop models (where q is the 
deformation parameter in the underlying quantum group symmetry) the periodicities 
of the boundary-related quantities e^^" and e^^'' is invariably twice that of e-^*" . On the 
other hand, for Ising models (where q enters the elliptic parameterisation of the 
coupling constant) the periodicities of all three quantities is the same. Another 
observation is than whenever the "same" model (Ising or Potts) is solvable on two 
different lattices, the periodicities are unchanged. It would be interesting to shed 
further light on these observations, for instance from the perspective of the corner 
transfer matrix and/or conformally invariant boundary states. 
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Table 1. Periodicities of the exponents appearing in the factorized form of bulk, 



surface and 


corner free enerj 


^ies 


for different loop models. 




Model 


Lattice 


e/b 


periodicity e-^= periodicity 


e-'^^ periodicity 


Potts selfdual 


Square 


4 


8 


8 


Potts selfdual 


Triangular 


4 


8 


8 


Potts antiferromagnet 


Square 


8 


16 


16 


Chromatic polynomial 


Triangular 


6 


12 


12 


FPL^ 


Square 


8 


16 


16 


Ising 


Square 


8 


8 


8 


Ising 


Triangular 


8 


8(?) 


8(?) 



3.3. Critical limits 

In this section we study the critical limit (s) g — t- 1~ (and g — t- — 1+ whenever the 
latter corresponds to a critical theory) of the product forms obtained for e^^'', e^^'' and 
e-^'^. Physically, one expects a finite limit for e^^'' and e^^", whereas e^^'' is expected 
to exhibit a divergent behaviour at the critical point, which will be related to a 
divergence of the characteristic dimensions of the system, such as its correlation 
length ^ ^. 

We first discuss a few general properties observed when studying the convergence 
of products having the structure ([2]), in the limit g — ^ 1^. We present these properties 
without formal proofs, but we demonstrate their applicability through a series of 
concrete calculations that we exemplify below. 

Let us first consider a product of the type ([2]) with purely periodical exponents 
ak (that is, with f3k = for all k), as it is the case for most of the expressions found 
in the course of this work. We can emphasize the periodicity property of such a 
product by rewriting it as 

oo 

JJ (1 _ qak-ia-yi _ ^ak-ia-y^ . . . _ ^aky^ ^ (g^^ 
k=l 

where a G N denotes the period. We have proved that this product has a finite hmit 
when g — )■ 1~ if and only if the two following conditions are satisfied 

71 + 72 + • • • + 7a = , (62) 
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7i(a-l)+72(a-2) + ... + 7,0 = 0. (63) 

If 71 + 72 + - • -+70 > (resp. < 0), the product is divergent (resp. has a zero hmit). If 
7i + 72 + • • • + 7a = 0, the product is divergent for 71 (a — 1) + 72 (a — 2) + ... + 7aO < 
and has a zero limit for 7i(a — 1) + 72(a — 2) + . . . + 7^0 > 0. 

Now turn to the general form ([2]), where = /S^fc + 7^. We have not found 
any solid analytical arguments enabling the classification of the g — ^ 1^ limit or 
the asymptotical behaviour of such products. However, the following property, yet 
unexplained, will turn out to be useful in the description of the Ising model corner 
free energy divergence: it was observed that products of the type 

°° /I „ak—m\ak—m 

n U ~ g ) /g.N 

1']^ Qak—n\ak—n ^ ' 

fc=l ^ ^ ' 

have a finite limit when g — )■ 1^ if and only if m + n = a (not mentionning the trivial 
case m = n, for which the product equals 1 for any q). 

3.3.1. Finite limits of hulk and surface free energies. To evaluate the limits of the 
e-^"" and e-^'' products, each factor (1 — g')"' can be replaced by (Z(l — g))"'. Using the 
identity 



k=l L ^ / 



k 



(65) 



where T{z) is the Euler gamma function, one is led to the following results. 



For the selfdual Potts model on the square lattice: 

lim e^-^ = 18—!^^^ ^ , (66) 

r(-i)^ 
r i'^Y 

lim e^= = 8 — . (67) 

r(-|)^ 

For the critical antiferromagnetic Potts model on the square lattice: 

lim e-^'' = , (68) 

g-s>l- 

3r(^r(|))- ^''^ 



Corner free energies and boundary effects 



29 



The g — i- 1^ (that is, Q ^ 0~) hmit of the bulk free energy is trivial, since in 
this case Z = both in finite-size and in the thermodynamic limit. However, 
dividing the finite-lattice partition functions Zmn by Q*^^ produces a model 
which in the Q ^ limit has a non-trivial combinatorial interpretation [3T| [32]. 
It amounts to counting so-called forests, which are collections of spanning trees 
(i.e., Fortuin-Kasteleyn clusters without cycles), on the square lattice, where 
each component tree carries the fugacity w = —4. The corresponding partition 
function is given by 

„ ^ e/b r(^)^ 

(^rorests)^/(^^")=lim^ = i)fi^. (70) 

g^i- Q r(|) 

Not much more about this problem seems to be computable from the expression 
for e-^^, since its derivatives with respect to Q involve clusters with cycles. In 
particular, it is unfortunately not possible to compute the average number of 
trees by this method. 

For the selfdual Potts model on the triangular lattice: 

-"^r(|)V(|f' 

lim e/= = 4 - 2^3. (72) 

For the chromatic polynomial on the triangular lattice [Ej: 

F (^^ 

lim e^'' = - 54 X 2^/^ — , (73) 
9^1- vr^F (3) 

lim = 8 . (74) 

vrF(i)^ 

For the FPL^ model: 

F (i)^ 

lim e^'' = — , (75) 
F (-^^ 

lim e/= = 144V27r ^ . (76) 



r(4) 



3.3.2. Corner free energy divergence. Using the previous method to compute the 
limit of e-^'^, one finds that the corresponding factorised form is either divergent or 
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has a zero limit. We shall see below that these two possibilities distinguish the sign 
of the central charge c of the CFT that emerges in the g — ?■ 1^ (or g — )■ —1^) limit. A 
more detailed study of the asymptotics of e-'^'^ will enable us to determine the precise 
value of c as well as the precise asymptotic divergence of the correlation length ^ in 
the critical limit. 

In this section we first expose the main tools needed in the asymptotic analysis, 
and we give precise results for each of the models for which we have been able to 
obtain an exact product formula for e-^'^ . In the following section 13.41 we confront 
these results with the CFT predictions 

To obtain the asymptotical behaviour of infinite products of the type ([2]), in the 
limit g — 7- one can use the properties of the Dedekind eta function, defined in the 
upper half of the complex plane by 

oo 

v{r) = n (1 - ' (77) 

k=l 

where q = e^^'^'^. Setting r = ix, and using the modular transformation identity 



Vi-r-') = V^vir) , (78) 
one can deduce the asymptotical behaviour of ri{ix) as g — t- 1^ (resp. x — ?■ 0+) from 
its Taylor expansion as g — (resp. x — ?■ +oo). Finally, we find the general formula 



oo 



nr 



1 —logo , 

— I ~ \l ^_ 66^. (79) 



k=l ^ 

In practice, we will use the related formula 



oo 



nr 



1 / 1-g , s 

r ~ \\oi ie6«(i-'?) , (80) 

_ - g°*-- q-^i- V 277 ' ^ ' 

from which the asymptotics of the different product forms obtained for e-^'^ can easily 
be obtained as follows. 

• For the selfdual Potts model on the square lattice one finds 



g/c 



k=\ 

oo 



9^1- 



nr^w)(nr3w) (n^g^^'^ 
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For the antiferromagnetic critical Potts model on the square lattice, we find by 
similar considerations 

For the selfdual Potts model on the triangular lattice the divergence of the corner 
free energy (l38l) . that is, the corner free energy of the three corners of angle | 
in a triangle-shaped lattice, is found to be 

fJ- ~ — e 6*^1-9 . (83) 

For the chromatic polynomial on the triangular lattice, the limit a; — ?■ — 1^ of 
the corner free energy (that is, g — )■ 1~ or Q — )■ 4"'') can be computed by turning 
the product fHOj) into a function of g, as indicated in f l44l) . This however leads to 
a large number of factors and lengthy calculations. We therefore indicate here 
an alternative way of evaluating the behaviour of the divergent contributions 
directly in terms of the variable x. In the product expression (HOl) for e^^'' the 
factors with an even power of x yield a finite contribution in the limit x — > — 1"*", 
whereas the diverging contribution is given by the factors with odd powers of 
X, namely of the type 1 — x^^*^"^^^'"'"^^ Series expanding the logarithm of these 
factors, permuting summations, and summing the geometric series over k we 
find 

°° °° °° ^(12fc-2p-l)?i 



iogn(i-x'='-*-') = -EE „ 

fc=l k=\ n=l 

^ ^-(2p+l)n ^ 

n=l 

- ^-(2p+l)n 
~ — > 



12n2(l 



TT 



X) 

2 



144(1 + x) 



(84) 



that is, in terms of q, 



71^ 



.ogn(l-(-r)— ■),^^^_^^. (85) 

independently of p. Summing up all the contributions, we find for the energies 
fc (l) associated with corners of angles | and ^ respectively the following 
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behaviours 



Mi 



7r\ Ttt^ 



3/ q^i- 4321 - g ' 

V 3 / q^i" 864 1 — q 

The hmit x — )■ 1~ (that is, Q — 0^) can be computed directly in terms of x, 
and we find using the same methods as for other models 

„ /7r\ Ttt^ 1 , , 

^=(3) „~--108— • 

. _Z^!^. (89) 
•'^ V 3 y 216 1 - x ^ ^ 

For the FPL^ model 

f 1 s^f./^^ 1\ 

qJ': _e 16 U-9 2) _ (QQ) 

q^l- 4 

For the square-lattice Ising model, the properties stated at the beginning of 
section 13.31 allow us to isolate the divergent contribution of the exponentiated 
corner free energy from finite factors. We thus find 



k=l 



g8fe-4)(l_g2fc-l)2 

~ ei6(i^ . (91) 

• For the triangular-lattice Ising model inscribed in an equilateral triangle, using 
similar methods, the divergence of e-^'^ in (l60l) is found to be 

e/c 0^ e 12(1-9) . (92) 
3.4- Relation to conformal field theory 

Our interest here is to understand these corner free energy divergences by the 
standards of conformal field theory (CFT). Cardy and Peschel showed in [3] that 
the presence of a corner of interior angle 7 along boundaries of typical size L in 
a two-dimensional conformally invariant (hence, critical) model of central charge c 
originates from a logarithmic correction to the free energy, namely 

AF= -^(l-(7r/7)2)lnL 
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^InL for 7 = 7r/2 , 

flnL for7 = 7r/3, (93) 
^InL for7 = 27r/3. 

Note the different sign compared to Cardy and Peschel's explicit formula, due to the 
fact that we chose to define the free energy as the logarithm of the partition function, 
without the conventional minus sign. 

We now wish to compare the results obtained above with this prediction. We 
emphasise the fact that the FLM formulae themselves forbid the appearance of In L 
like terms, since all the contributions /b, /s and /c are by construction independent 
of the system's dimensions. Indeed, at the critical point (? = 1, the series expansion 
in powers of q breaks down, since all the finite subgraphs so far neglected by the FLM 
formulae become non-negligible. This problem can be reformulated in other terms: 
if the M X N lattice is considered finite, Enting's formulae give an exact result up to 
order k ^ M + N = L. When approaching criticality (g — i- 1^), the cutoff B{k) has 
thus to be chosen bigger and bigger in order to reach a given precision in the series 
approximation, and so must therefore the size of the system. So to summarise, our 
use of the FLM formulae is valid at criticality only provided that the size L of the 
system diverges fast enough to include excitations of increasing size, that is, diverges 
as the correlation length. In this context eq. ( I93l) can thus be rewritten as 

AF = -^(l-(7r/7f)lne, (94) 

where ^(g) is some characteristic length defining the typical maximal size of 
excitations that enter the FLM formulae with a non-negligible statistical weight. 

3.4-1 ■ Selfdual Potts model. For the selfdual Potts model on the square lattice, we 
have c = 1 when approaching the critical region via the limit g — i- 1~. Recalling that 
the corner free energy divergence is the sum of four 7 = 7r/2 contributions of the 
type (pSi) . we thus have 

q^i- 2(1 - q) 

or more precisely, 

1 
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It is of course tempting (as suggested by our notation) to interpret ^ as 
a correlation length for the near-critical system. We shall now see that this 
interpretation is indeed correct. Based on Bethe Ansatz (BA) calculations, Wallon 
and Buffenoir [3l] have identified the correlation length ^ba from the ratio of the 
leading and next-leading eigenvalues of the XXZ spin chain hamiltonian. They obtain 
asymptotically 

^BA ~ -^e^. (97) 
8v2 

This coincides with ( 196|) up to a g-independent multiplicative prefactor (of numerical 
value ~ 1068). Crucially, the constant ^ in the exponential, governing the strength 
of the essential singularity, is precisely the same. The multiplicative prefactor was of 
course to be expected, since a system does not have just one correlation length, but 
several, which are all proportional and thus present the same critical divergence. We 
can therefore safely refer to ^ as a (and sometimes, by an abuse of language, even 
"the") correlation length. 

We now turn to the selfdual Potts model on the triangular lattice, for which 
eq. ( l83l) provides the asymptotics of the corner free energy for three corners of angle 
|. In conjunction with eqs. ( I93|) -( p4|) . and using that c = 1 is lattice-independent, 
we obtain 

^ 8e^'/4 ^ ' 

The fact that the universal, g-dependent part of the asymptotics precisely coincides 
with ( 19 6 p for the square lattice is a nice verification of universality, and of the angular 
dependence appearing in the Cardy-Peschel formula (1931) . 

Note also that ratio of the correlation lengths on the two lattices tends 
asymptotically to a constant, 

lim g ^ 1" 1^ = 128(3^3 - 5)^ ~ 0.966031. (99) 

It would be interesting to confront this prediction with numerical results, such as 
Monte-Carlo simulations. 

In section H] below we shall investigate how the effective central charge — 
extracted from the divergence of the corner free energy using (|93|) — is affected by 
changing the boundary conditions at the position of the corner. This change can be 
intrepreted within CFT as the insertion of a boundary condition changing operator 
in the corner. 
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3.4-2. Antiferromagnetic Potts model. For the antiferromagnetic Potts model on 
the square lattice, the limit g — t- 1^ (or Q — )■ 0^) corresponds to a critical theory 
with c = — 1 lini HH [15]. The comparison between f l82|) and fl94|) gives the following 
prediction for the behaviour of the correlation length 

^(g) ~ e'^'^^"^). (100) 



As already mentioned in section [2^ numerical evidence has been given [TT] that 
the same c = — 1 CFT arises as the continuum limit of the triangular-lattice selfdual 
Potts model, in the limit where t — )• — 1"*", i.e., {Q,v) — )■ (0, —3). The corresponding 
correlation length can be computed from fl38l) . which is shown to diverge for t — )■ — 1+ 
hke 

e^'^ ~ e'^iiTt-^) . (101) 

Comparison with flM|) yields for the associated correlation length 

^(t) ~ e'^(TTi-5). (102) 

In order to compare the correlation lengths fllOOl) and (11021) . we need to express 
both in terms of the same parameter, most naturally Q. We find from the respective 
definitions of q and t that 



2(1 -g) = ^(l-^^+OiQ)) 



(103) 



for the antiferromagnetic square-lattice Potts model, and 

(104) 



3(1 + t) = ^(l-^^+0{Q)) 



for the triangular-lattice selfdual Potts model. Both results can thus be rewritten in 
the common form 

^(t) ~ . (105) 

This provides a verification of the universality between the two models, which goes 
beyond the numerical evidence of [llj. Note that not only do the two correlation 
lengths given above have the same critical divergence, but the prefactors are also 
identical (and equal to unity). 
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3.4-3. Chromatic polynomial. For the chromatic polynomial we can perturbatively 
access two different critical theories. In the limit g — t- 1^ (or Q — )■ 4"*") we have a 
CFT with = 2, whereas the limit g — t- z x 1~ (or Q — 0^, or a; — ?■ 1^) one finds 
c = —2 (see [22j and references therein). In both cases we must pay attention to the 
fact that the corner angles in have to be adapted to the triangular geometry. 
For the g — ?■ 1~ limit we find 

96(1^ 

for a corner of angle |, and 

for a corner of angle Meanwhile, for the g — )■ i x 1~ limit we find as a function 
of X 

for a corner of angle |, and 

15(1^ • (™' 

for a corner of angle 

It is obviously worrying that for both limits the results for the divergence of 
In^(g) depend on the corner angle, i.e., are mutually inconsistent. On the other 
hand, we have seen that both the results for the selfdual Potts model and those for 
the antiferromagnetic Potts model nicely confirm universality between the square and 
triangular lattices, once the angular dependence of (!93|) has been taken into account. 
(We shall also obtain a similar agreement for the Ising model in section [3.4.5l below.) 

One possible explanation for the discrepancy is that the corner angles seen in 
the continuum limit do not equal those on the lattice (i.e., | and ^). This might 
be the case for geometrically highly constrained (frustrated) problems, such as the 
4-colourings under consideration. 

An example of a situation where such a scenario is known to hold true is provided 
by dimer coverings of the so-called Aztec diamond lattice, which is just a diagonally 
oriented square piece of the square lattice. It has been proved rigorously [35] that 
there are Z = 2^(^+^)/^ dimer coverings of an Aztec diamond of size A^. It follows 
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that /c = exactly for any A^. Meanwhile, it is well-known that dimer coverings 
of the square lattice are in bijection with configurations of a scalar height, whose 
continuum limit is a free Gaussian field. Therefore the central charge is c = 1. This 
would seem at odds with the result fc = 0, since the Cardy-Peschel result ( l93l) then 
predicts that fc depends non-trivially on A^. The resolution of the apparent paradox 
is that, due to the boundary conditions, the dynamics of dimers close to the corners 
of the lattice is frozen. In the thermodynamic limit, the region where dimers are free 
to move becomes — by the so-called Arctic circle theorem [36] — a circle inscribed in 
the Aztec diamond. In this sense there are no comers in the continuum limit, and 
the paradox is resolved. 

We believe that it would be interesting to investigate whether the 4-colouring 
problem has frozen dynamics close to the corners in the thermodynamic limit. 

3.4.4. FPL'^ model. The FPL^ model has central charge c = 3 when g -> 1~ [26l, [27]. 
We thus find from that 

lne(9) ~ jn^v (110) 
q^i- 4(1 - q) 

3.4.5. Ising model. 

• For the Ising model on the square lattice, of central charge c = | at criticality. 



we find 



Ine(g) ~ 7^. r- (111) 

2(1 - q) 

This result can be compared to the analytical expression given in [18] for the 
Ising model correlation length in the critical limit, namely 

r'~-lnfc, (112) 
where k is defined as 

k = {smh2K)-^. (113) 
Rewriting this as a function of the temperature T, we have 
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from which one infers the value z/ = 1 of the correlation length critical exponent. 
In order to retrieve this well-known result from f illip , we need to re-express the 
latter in terms of T. 

It is easily checked that below the critical temperature one has 

Xc — x oc Tc — T . (115) 

We thus need to compute the asymptotical behaviour of x in terms of q when 
q — )■ 1~. This is done in terms of elliptic functions of modulus k from eq. (3) of 
[29] . namely (with the usual notations) 

iK' 

= -ik^/^n— , (116) 

where q = exp (-^). The limit g — 7- 1 is obtained for /c — )■ 1 , yielding for K 
and K' the following expansions 

K ~ -lln(l-A;2) _ + In4 + O (l - A;^) , (117) 



^' - i (^1 + ^ + ^ (1 - k-f) ■ (118) 

Series expanding eq. f lll6p around fc^ = m — )■ 1^ yields 

= -1 + 72 - 8 (2 - v^) (1 - + C ((1 - ky) , (119) 
that is, to leading order, 

Xc-x oc Xc^-x^ oc -(1-A;^). (120) 

From the above definition of q in terms of K and K' we furthermore have 

1-q ~ — — , (121) 

^ k^i- 41n(l-A;)' ^ ^ 

and thus, 

Tc-T ~ (l-k)'^ ~ e"^^. (122) 

We can thus rewrite the divergence (11110 of the correlation length in terms of 
Tc — T, which yields, as expected, 

e (X p^T^- (123) 
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For the triangular-lattice Ising model, in the case of a lattice inscribed in an 
equilateral triangle, we find from ( l92l) that 

In the same way as for the square lattice, we seek an expression of the parameter 
X defined in (15^ in terms of elliptic functions. One possibility is given by 

a;2 = -ikkn (^-^^ , (125) 

where, once again, q = exp (— Using ( IllSp . we find that around the critical 
point Xc"^ = 

that is, as in the square-lattice case, 

Tc-T ~ (l-kf- (127) 

We thus have once again 

e oc ,^ , (128) 



T-s>Tc 



\T-Tr 



which gives in comparison with f ll23p a verification of universality between 
square and triangular-lattice Ising models. 



4. Particular boundary conditions 

We now consider the effect of particular (non-free) boundary conditions on the 
previous models. In section 14. 1^ we derive FLM formulae for the case of lattices 
where one or more sides are endowed with particular boundary conditions. These 
formulae allow us to compute the surface free energy with the particular boundary 
conditions, or the corner free energy for a corner between two sides with free- 
particular or particular-particular boundary conditions. The resulting expressions 
will be compared with those found in section [3] for the case of free boundaries 
boundary conditions. 

As will be discussed in further detail in section 14. 3[ the Ising model presents 
some special difficulties, whose resolution we leave for future work. Therefore, we 
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will focus for the remainder of this paper on the selfdual Potts model on the square 
lattice|§| The particular boundary conditions that we shall be interested in are those 
introduced by Jacobsen and Saleur [7]. We refer to them as JS boundary conditions 
for short, and define them precisely in section S31 

Our results for JS boundary conditions are presented in section IT5l The critical 
limit g — )■ 1~ is studied in section 14.61 where we also make contact with the CFT 
results of [7j . In particular, we shall see how physical observables such as correlation 
length and conformal properties are modified by the change of boundary conditions. 

4.1. FLM for particular boundary conditions 

In the case of lattices with particular boundary conditions on one or several sides, 
the finite-lattice expansions given in section [3] must be modified to take into account 
new types of contributions, where the sublattice touches one or more sides of 
the [M, N] lattice. 

In [30j, Enting has shown how to derive numerically FLM formulae in the 
case of fixed boundary conditions for the 3-state Potts model. In this section, we 
demonstrate how particular boundary conditions can be dealt with analytically in 
a framework that does not depend on the explicit choice of model or boundary 
conditions. In particular, we shall derive explicit formulae for the different 
configurations corresponding to one, two, three our four sides with particular 
boundary conditions, regardless of what these conditions might be. 

By convention we choose to denote the sublattice contributions as follows: 

• fm,n for an m X n lattice with free boundary conditions; 

• fm!n (resp. for a particular boundary condition on one vertical (resp. 
horizontal) side; 

• /mil particular boundary conditions on two adjacent sides (i.e., one vertical 
and one horizontal); 

• fl^^ (resp. f^n) particular boundary conditions on two opposite vertical 
(resp. horizontal) sides; 

§ We could equallly well have studied any of the other Potts or fully-packed loop models with 
particular boundary conditions. But technically it is easier to tackle the square-lattice Potts model, 
since in this case 1) the coefficients of the relevant product formulae present the smallest periodicity, 
and 2) our numerical methods for computing finite-lattice free energies are more powerful. 
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• (resp. f^n) particular boundary conditions on three adjacent sides; 

• f^n particular boundary conditions on all four sides. 

Similar notations apply to the / contributions. 

Just as in the case of free boundary conditions, we start in the case of a bounded 
M X N infinite lattice by writing the decomposition 

fit= E E (#[M./3],[M,iV,„])4?\ (129) 
i<M,j<N 13 

where #[j,j,/3],[M,7v,a] denotes the number of ways that an i x j lattice with 
boundaries can fit into the considered M x N lattice with (a) boundaries. The 
self-consistent equations relating the finite-lattice contributions f^j'' and fm)n are to 
be written in the same fashion (see below). 

4.1.1. Particular boundary condition on one side. For one particular boundary 
condition on (say) the left side, f ll29p takes the form 

fM,N= E (M-^)(iV+l-j)/M 

i<M-l,j<N 

+ iN+l-j)fy, (130) 

i<M,j<N 

where the first (resp. second) term corresponds to contributions where the [i,j] 
rectangle does not (resp. does) touch the left side of the [M, A^] rectangle. 

The fij were computed in section |31 In the same way, the flj are determined 
self- consistently from (11301) — with the M x N lattice being replaced by an m x n 
finite lattice — and are found to be 

X]= E V{n,i){Sm,^- Sm,i-l){fl;^^^- fm-l,n) ■ (131) 
m<i,n<j 

Introducing the cutoff set B{k), the free energy of the bounded infinite lattice is thus 
approximated as 

fM,N- E fm,n{M-t){N + l~jUm,t)v{n,j) (132) 

lm,n],[i,j]<B(k) 

+ E - fm-l,n){N + 1 - j){5m,t " 5m,i-l)vinJ) ■ 

[m,n],[i,j]<B{k) 
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After performing the summation over the sublattices and having isolated the 
free boundary contribution from the corrections related to the particular boundary, 
we have 

& = fM,N + iV5/P + 5/P , (133) 



where 



and 



E (fi-/m-l,n)x 

[m,n\<B{k) 

n-1 + 5m,A:-n-2) — /b (134) 



^/c ^ E (/U-/r«-l,n)(m-fc + l)x 

[m.,n\<B(k) 

)-fs. (135) 

Notice that the bulk energy is left unchanged by the introduction of the boundary, 
and that the surface correction is logically proportional to the length of the 
corresponding side. Similar remarks will hold with two or more particular boundaries. 

4.1.2. Particular boundary condition on two adjacent sides. In the same way, we 
have for two adjacent particular boundaries 

fM,N= E (M-z)(iv-j)A,+ E 

i<M-l,j<N-l i<M,j<N 

+ E (N~j)fV+ E (M-t)f^, (136) 

i<M,j<N~l i<M-l,j<N 

with the fij and j already detemined above. Self-consistency of eq. (11361) for finite 
lattices requires that 

^ — 5^ {Sm,i — Sm^i-l){5nJ — Sn,j-l) X 
m<i,n<j 

(/m,n "I" fm~l,n~l ~ fm,n-l " /m-l,n) • (137) 



Finally, 



& = fM,N + (M + N)SfP + , (138) 
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where 



^fc (/m,n fm-l,n-l fm,n-l fm-1, 

lm,n]<B(k) 



{Sm,k-n - Sm,k-n-l) + - 25/^^ - /b . (139) 

4.1.3. Particular boundary condition on two opposite sides. For particular 
boundary conditions on two opposite sides (say, the right and left sides), similar 
calculations lead to 

& = fM,N + 2N6fP + , (140) 



where 



(5/P = 25/P. (141) 



4.1.4- Particular boundary condition on three sides. For particular boundary 
conditions on three (say, the right, lower and left sides), 

& = fM,N + {2N + M)SfP + 6f^ , (142) 



where 



5/P = - 5/P . (143) 



4.1.5. Particular boundary condition on all four sides. Finally, for particular 
boundary conditions on all four sides, 

& = fM,N + (2iV + 2M)6fP + dfP , (144) 



where 



5fP^A5fP-A5fP. (145) 



4.2. General observations 

The remarks made in the case of can be extended to all the boundary 

configurations. In general, the above calculations show that: 

• The bulk free energy is not affected by any change of boundary conditions. 
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• The surface energy can be written as a sum over independent sides, the energies 
associated with each type of side being distributed according to 

fi = /s/2 

/; = /! + 5/1, with 6 fl^ 6 fP. (146) 

• The corner energy can also be written as a sum over independent corners, the 
energies associated with each type of corner being distributed according to 

f\r = /P/4 

fp = fir + l^fP (147) 

fP = f\r + 6fP-6fP. 

The interaction energy between two adjacent sides of the lattice is taken into 
account by the corner energy. However, one does not observe any corner-corner 
interaction term, which is related to the assumption implicitly made in the FLM 
formalism that the MxN lattice is bigger than any cutoff set B{k), that is, effectively 
infinite. 

The absence of corner-corner interactions is radically different from the outcome 
of CFT calculations, where the partition function on a large rectangle depends non- 
trivially on its aspect ratio through the so-called modular parameter. This implies 
in particular (but not only) that the change of boundary conditions at the corners 
interact, and in fact define a correlation function of boundary condition changing 
operators. In view of this fundamental difference, it is all the more remarkable that 
the asymptotics of the corner free energies reported here link up nicely with CFT 
predictions 



4-3. Case of the Ising model 

As previously announced, the case of the Ising model with fixed boundary conditions 
is not correctly described by the above calculations. 

To get a feeling of what goes wrong, we first consider the example "H — ff", 
where the two adjacent sides of the MxN rectangle support respectively + and — 
fixed boundary conditions, and the remaining two sides are free. This implies that a 
domain wall will originate from the H — corner and terminate on any of the two free 
sides. In particular, its length will be at least min(M, A^). This situation is at odds 
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with the perturbative principle illustrated in Fig. |5l where the excitations are with 
respect to a trivial ground state. 

So contrary to what was done in the case of loop models, finite-sublattice 
excitations in the Ising model with one or more sides supporting fixed boundary 
conditions cannot be thought of as localised perturbations independent of the 
surrounding ground state. Instead, these excitations interact with the surrounding 
spins, therefore affecting their configuration of minimal energy. We leave the 
resolution of this problem for future work. 

In the following, we shall therefore focus on boundary conditions in loop models, 
where excitation can still be described as local perturbations of the ground state in 
Fig. [5l This leads us naturally to consider the JS boundary conditions [7] for the 
selfdual Potts model on the square lattice. 

4 ■4- JS boundary conditions 

Within rational CFT the number of possible conformal boundary conditions is known 
to be equal to the number of primary operators. This result does however not apply 
to the Q-state Potts model for generic values of Q. Instead one expects the existence 
of infinitely many distinct conformal boundary conditions. 

One infinite family of such conformally invariant boundary conditions was given 
in [7]. Let us parameterise the bulk loop weight n by 



so that p G M reduces to the usual minimal model index whenever p G N. The JS 
boundary conditions then amount to assigning a different weight ni to each loop 
touching at least once any one of the particular edges. Following [7] we parameterise 
ni by 



with < r < p + 1. 

In the following we shall only be concerned with the case of r G N. The case 
r = 1 corresponds to Ui = n, or Dirichlet boundary conditions, and the corrections 
to surface and corner free energies should therefore be exactly zero. When r = p/2 




(148) 




(149) 
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we have Ui = 1, which corresponds to Neumann boundary conditions. Note also 
the case r = which corresponds to rii — ?■ oo. Each spin on a particular boundary 
is then enclosed by a distinct boundary loop, a problem which seems likely to have 
physical applications. 

The parameterisations fll48l) - fll49p need to be linked up with the fact that the 
FLM method permits us to approach the conformal limit n — t- 2+ from the side n > 2 
(i.e., Q > 4). We therefore set 7 = and let 7 = 17, with 7 G IR+. Thus 



In the conformal limit 7 — )■ 0''' we thus have Ui ~ meaning in particular that 
ui > 1 for any finite r G N, whereas the Neumann boundary condition ni = 1 is 
formally obtained in the limit r — )■ 00. These remarks will turn out important for 
the following discussion. 

4.5. Results for JS boundary conditions 

In this section we report our results for the corrections to the surface and corner free 
energies for the selfdual square-lattice Potts model with JS boundary conditions. 
We have performed the explicit FLM calculations for r = 0, 1, . . . , 9 and from those 
we have been able to conjecture product formulae that we believe are valid for any 
r G N. 

4.5.1. Bulk free energy. The bulk free energy is invariably found to be unchanged 
upon taking r 7^ 1 in the JS boundary condition. This was of course to be expected 
on physical grounds, and indeed follows analytically from the results of section 14.21 

4.5.2. Surface free energy. The factorised form of the correction to the surface free 



energy (i.e., 5f^ = fs ~ fl with = f^^ and fl = f^^"^^) for the selfdual Potts model 



n = 2 cosh 7 , 

sinh((r + 1)7) 



(150) 



on the square lattice reads, as a function of r: 



• for r = 0, 




(151) 
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for r = 1, 

e-^^i = 1 , 

for r an even positive integer, r = 2p. 

p-i 

m 



1 -g4^+i)(l -g^ 



g4fc-3)(l _g4fc+4) 

for r an odd positive integer different from 1, r = 2p + 1, 

.■5/1 



e 



1 (l_g4fc-l)(l_^4fc+4) • 

When r — 7- oo, f|153p and f ll54p give the same hmit, namely 



4-5.3. Corner free energy. In the same way we find the corrections to the 
free energy when one side supports the JS boundary condition: 

• for r an even positive integer, r = 2p, 



n X — n'^k+i n (I 



_ q2k+l ii (I_g8fc-l)(l_g8fc+l) 

oo 

ll(^-q''^'){l-q''-'), 

k=l 

for r an odd positive integer, r = 2p + 1, 



- — q 

k=2p+l ^ k=l 

which includes the special case e ■''^ = 1 for r = 1. 
When r — > oo, the common limit of f ll56l) -( fT57|) is 

k=l 

Similarly, when JS boundary conditions are imposed on two adjacent sides: 
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• for r an even positive integer, r = 2p, 



□ 1 - g4p+2 

e ■'^'^ = 2—— X 



1 -g" 



^^^^^ g4fc-l)2(l _ g4fe-3)2(l _ g8fe-6)(l _ ^8/=-2) 

«^ /I „8fc-5N2/i „8fc-3N2 

(l_g8A:-2)(l_g8fc+2) ' ^"^^^^ 



for r an odd positive integer, r = 2p + 1, 



5fP _ 1 - q^"^^ fr (l-g4fc-l)2(l_g4fc+l)2 



(1 - g4fc+2)2 ' (160) 



fc=l 



which includes the special case e ■''^ = 1 for r = 1. 
When r — > oo, the common limit of fll59l) - fll60p is 

k=l 



4-6. Critical limit 

In the critical limit, g — )■ 1^, we can now extract the finite limits of the surface free 
energy corrections, and the asymptotic divergent behaviour of the corner free energy 
corrections. 

4-6.1. Finite limits of surface free energy corrections. The g — t- 1 limit of e is 
found to be: 

• for r = 0, 

linie^^! = -^^, (162) 
vrr(i)^' 
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• for r an even positive integer, r = 2p, 

lime^^! = 8V2^ r(|+p)r(i+ri 



r(i)^r(i+p)r(| + p) 

for r an odd positive integer, r = 2p + 1, 

^sf\ ^ 2 r(|)r(|+p)r(|+p) 
0Fr(f)r(f +p)r(2 + p) ' 

which includes the special case liniq^i- e = 1 for r=l. 
When r is taken to infinity. 



(163) 



(164) 



(165) 



4-6.2. Divergence of the corner free energy corrections. More interesting is the effect 
of JS boundary conditions on the critical divergence of the corner free energy, which 
gives access to the effect of these particular boundary conditions on the effective 
central charge. 

In the case of a corner between a free edge and one edge with JS boundary 
conditions, the correction to e-^'^ is found to have a finite limit 



for any finite value of r (whether even or odd). This finite limit just adds a constant 
term to the corner free energy, with no effect on the diverging part. 

On the contrary, when r is taken to infinity, we find from fll58p that 



e-^^P ~ Jl-^e"^^^"^^ (167) 

so that the diverging part is modified. We shall discuss these findings more fully 
below. 

In the case of a corner between two edges with JS boundary conditions we find 
for finite r 

lin, e^^P = ^2-^V2 + v^ vr3/^r (If^f 

,^1- 2 V r(^)'r(^)'r(^)' 
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On the contrary, when r is taken to infinity, we find from f ll6ip that 

e-s/c ^ 1 e-24[—-2) . (169) 

q^i- 4 ^ ^ 

Note that, apart from the prefactor, the diverging part is identical to that of ( I167p . 

In conclusion, the diverging parts of the corner free energies associated with 
each possible kind of corner are 



71^ 



,^1- 32(1 -g) 



~ — ; for r — oo . 

q^i- 32(1 -g) 

4-7. Relation to conformal field theory. 

We can now compare the divergence of the corner free energies with CFT results for 
the JS boundary conditions |7]. 

It was found in |7] that the operator that changes the boundary conditions from 
free to JS with parameter r has the conformal weight 

^2 _ I 

h = hrr = -r-, TT. (171) 

In the limit n — 2^ (i.e., "p — )■ oo) any fixed boundary loop weight n\> \ corresponds 
to a value of r that grows slower than p. In other words, /i — )■ as p — )■ oo. On the 
other hand, n\ = \ corresponds to r = p/2, in which case we obtain the finite limit 
— 7- ^ as p — 7- oo. So summarising, h{n\) tends to a step function: 

{\ forni < 1 , 
-i. for ni = 1 , (172) 
ID 
for rii > 1 . 

On the other hand, the insertion of a boundary condition changing operator in 
a corner will change the effective central charge according to [37] 

Cgg = c-32/i. (173) 

So according to (!93|) we expect a change in the divergence of the corner free energy 
if and only if /i / 0. 
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Recall meanwhile that ni ~ from fllSOp . So whenever r is finite we have 
TT-i > 1, and so by fll72p the divergent part of the corner free energy should remain 
unchanged. This is precisely what we have found to be the case in section I4.6.2I 

On the other hand, in the r — )■ oo limit we have ui = 1, and so from f ll72p -( |T73l) 
the divergent parts of and f\r- should be different. While such a difference is 
indeed apparent in eq. (I170p . the actual values do not quite work out as expected: 
we have c = 1, and from f ll70p we find Ces = |, so that h = This is at odds with 
h = jq given in fll72p . 

This discrepancy is maybe not completely surprising. Indeed, two very different 
double limits are at work in the conformal case (p — )■ oo and rii — )■ 1) and in the 
asymptotic analysis of the corner free energy (g — 1~ and r — > oo). Meanwhile, the 
critical exponent tends to a step function fll72p . so the only disagreement concerns 
the value of h right at the step. It is conceivable that a non-commutativity of limits 
misses this value (by a factor of 3). 

It remains to discuss the case of /|~. Since in this case there is no insertion of 
a boundary condition changing operator in the corner, the CFT prediction fll73p is 

that /|~ should have the same divergence as /|^. This is precisely what we have 
found in fll70p . On the other hand, it is clear that the two types of corners are 
different, and this should be reflected by a finite (non-diverging) difference in the 
two corner free energies. Once again, this is exactly what we found in section 14.6.21 

5. Conclusion 

In this paper we have presented what is — to our knowledge — the first systematic 
study of corner free energies from the perspective of exactly solvable models. 
Combining the FLM method with exact enumeration results, we have obtained 
exact (albeit conjectured) product formulae for the corner free energy e-^'^ for several 
integrable cases of two-dimensional Potts and Ising models, as well as for the FPL^ 
loop model. 

We have obtained the asymptotic expansions of these expressions for e-^"" near 
several conformally invariant critical points. This has permitted us to identify the 
asymptotic divergence of the correlation length, in agreement with Bethe Ansatz 
results whenever the latter are available, and provided new results in other cases. 
More importantly, the comparison between results for the square and triangular 
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lattices gave agreement with ideas of universality and has enabled us to verify the 
angular dependence of the Cardy-Peschel formula (l93l) . Such agreement was found 
in particular for the ferromagnetic and antiferromagnetic transitions of the Potts 
model, and for the Ising model. 

In some cases we have also provided new results for the surface free energy e•^^ 
In the first part of the paper we were concerned with free boundary conditions. 
But in section H] we have shown how to generalise the FLM formalism to the case 
of special boundary conditions. We have used this to study in detail the so-called 
JS boundary conditions [7], taking the selfdual square-lattice Potts model as an 
example. In particular we found some agreement with CFT predictions for the case 
where the corner contains a boundary condition changing operator. 

We leave several issues for future work. On the side of exactly solvable models, it 
would be interesting to establish if the corner free energy can be obtained exactly from 
the corner transfer matrix. If this is possible, one could hope to prove our formulae 
for e'^''. On the conformal side, we believe that our results can be interpreted in terms 
of boundary states, generalising the ideas of |37] . 
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Appendix A. Factorised form of the bulk free energies from analytical 
results 

In this appendix, we describe how certain factorised expressions found for the bulk 
free energy of Potts and Ising models can be recovered from existing exact results. 
The first step is to rewrite into a slightly different form the logarithm of generic 
factorised expressions of the form 

oo 

■Q(l _g/3fc+7)Mfc+^'. (A.l) 
k=l 

Using a series expansion, we have 

oo oo 

log - =Y^(^k + v) log (1 - g^'=+^) 

k=\ k=l 
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oo 



^ oo 

- ^ (/iA; + u)q 



^ n 

n=l k=l 



E 



n=l 



n 



•,/3r^^2 



(A.2) 



where between the last two hnes we have used the exponential polynomial summation 
formula to perform the sum over k. From this formula, we can identify some of the 
bulk free energies obtained in the main text with existing exact results. 



Appendix A.l. Selfdual square-lattice Potts model 



Applying the above identity (IA.2P to the product form fl271) obtained for the bulk 
free energy of the selfdual square-lattice Potts model yields 



oo ^ oo ^ 

/. = 21ng + ln(l + g^) + 25^|-45^-- 



,3fc 



k=l k=l 
^ -,k 1 ^2k 



+ q 



2k 



21ng + ln(l + g2) + 2 V^- 



k=l 



q 1 — q 



2k ■ 



(A.3) 



In [18], the dimensionless bulk free energy (that is, the opposite of our /b) of 
the square- lattice Potts model is shown to be ( l28l) . which we reproduce here for 



convenience, 



-lnQ-2 
2 ^ 



oo _ 



-kX 



k=l 



sinh 2k P 
cosh kX 



where P and A are defined by 

Q^/^ = 2 cosh A 
V sinh P 



QV2 sinh(/3 - A) ' 
On the selfdual curve, v = a/Q, that is, /3 = | with —A = Ing. We thus have 



sinh2fc/3 
cosh kX 



2k 



tanh kX 



1-g 



and thus 



V' = -ln(g+-)-2 



1 q 1 

k=l 



q 1 — q 



2k 



2k 



(A.4) 

(A.5) 
(A.6) 

(A.7) 
(A.8) 



that is ip = —fh as claimed. 
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Appendix A. 2. Antiferromagnetic square-lattice Potts model 



For the critical antiferromagnetic Potts model, our aim is to approach the critical 
point (Q, v) = (0, 0) from the left side of the v = —2 + a/4 — Q curve, that is with 
Q < 0, using the parameterisation 



Q 



(A.9) 



Assuming that /b has no singularity at Q = oo, this region can be mapped to the 
Q > 4 region by setting 

q = ip, (A. 10) 

and taking p real (more precisely, p G [0, 1[). Since |g| < 1 all along the problem, 
the FLM expansions are kept convergent and should remain valid. In terms of this 
parameterisation, we have on the v = —2 + curve 



V = — 2 + i{p — p~^ 
We now can use Baxter's formula in 



namely 



where 



K 



- In Q - 2 

2 ^ 



u 



-k\ 



k=l 



sinh 2ku 
cosh k\ 



(A.ll) 



(A.12) 



— fh , 

2 cosh A , 
sinh(A + u) 



(A.13) 



sinh(A — u) 

In terms of p, we find from (lA.lip and (]A.14p the following expression for u, 



u 



Sivr 1 , 

In » 

4 2^ 



Sivr A 
-+2 



(A.14) 



which indeed satisfies the conditions < < A and < < tt. Eq. ( ]A.12p can 
then be rewritten 



- In g - 2 
2 ^ 



Sivr 1 ^ i-ipf 1 - (-l)V 
In » + > — — — 



- ln(l +p^) + 21nj9 



Sivr 



k=l 

oo 

k=l 



+ p 



l+P 



2k 
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ln(l - q^) + Ing^ 



2 ^ k 1 + g2fc 



2fe 



Writing the sum on the right-hand side as 



y- 

^A: 1 + g 



2fc 



fc=l 
oo 



k 1 - g4fc 
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(A.15) 

(A.16) 



A; 1 - g4fc 2^ A; 1 

fe=l ^ k=l 



q^^ q^ 



T—- 
^ k 1 



fc=i 



we see from flA.2p that it can be put in the form 



In 



;i-^)nr 



fe=l 



-q- 



(A.17) 



(A.18) 



Inserting this back into ( 1A.15I1 . e~'^ is finally shown to be equal to the factorised 
expression ( IHT]) of ^ up to a — i phase factor which is presumably due to the 
crossing of some branch cuts during our analytic continuation. 

Appendix A. 3. Self dual triangular-lattice Potts model 

For the selfdual Potts model on a triangular lattice, the free energy in the 
thermodynamic limit is found in [18] to be given by ( 135|) . which we reproduce here 
for convenience, 



V' = -^lnQ-3 



k=l 



sinh 2k P 
cosh kX 



(A.19) 



with the same notations as in section 13.2. 1[ In terms of g, we have on the selfdual 
hue (ED 



e-^ = g, 

such that eq. (135|) reads 
1 



V' = --ing-3 



--\nq + Y^ 



^ q2k q-2n/3 _ q2n/3 



k=l 



l + q 



2k 



(A.20) 



(A.21) 
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Applying (]A.2p to the factorised expression (1341) for e-^*" , we check that it agrees with 
the above result, that is, —iIj = f^. 

Appendix A. 4- Ising model on the square lattice 

Applying eq. ( 1A.2I) to the expression ( H9l) found for e^^'^ in the case of the square- lattice 
Ising model yields 



/b 



1 °° 1 
- In (7 — > — 

2 ^ 



+ 



n=l 

g^"(2 - 2gS") 



8n^2 



[1-q 



8n^2 



8n\ 



1 - q 



8n^2 



8n^2 



8nN 



.1 - g 



g^"(5 + 3g' 
(1 - g8")2 



^2ra 



^6n 



1 _ g8n ]^ _ 



8n 



--lng + ^ 



n=l 



l + g"(-l + g2 + g 



2n\ 



nl + q 



An 



2n\2 



(A.22) 



In the same fashion, Baxter's analytical expression (]53l) can be transformed as follows 

,n\2/i „2nN2 



^ nil 



n=l 



^Ing + lnH}^ 



+ g2ri)(l - 



1) ^ n(l 



^n\2/ 



^2n\2 



IJL (l_g8fc-7)(l_^8fc-l) ' n(l+g2n)(l_g8n) 



1 

2^"^ + i-n 



n=l 



l + g"(-l + g2 + g 



2n^ 



1 + g 



4n 



2ri^2 



(A.23) 



which establishes the equivalence of the two expressions. 



Appendix A. 5. Ising model on the triangular lattice 

Applying (1A.2P to the expression (13^ found for e-^*^ in the case of the triangular- lattice 
Ising model yields 

Q °° Sn 

2 ^ n(l - g8") L ^ ^ ^ J 



^ nil 



-,8n 



n=l 



q 



[(6-3(l-g8")(g-^^"/3_^-i0n/3)^ 



8n\2 



+ [6(g"^"/^ + g~2n/3 _^ ^8n/3 _ ^-8n/3 _ ^2n/3 _ g4n/3^)j 
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oo 



3g^" 1 - q 



,2n 



n(l - gSn) ij^q. 



,2n ' 



(A.24) 



n=l 



which is seen to be equivalent to the analytical expression f l57|) derived by Baxter. 
References 

[1] H. Dosch, Critical Phenomena at Surfaces and Interfaces, Tracts in Modern Physics 126 
(Springer, Berlin, 1992). 

[2] J.G. Brankov, D.M. Dantchev and N.S. Tonchev, Theory of critical phenomena in finite size 

systems, (World Scientific, Singapore, 2000). 
[3] E.K. Sklyanin, J. Phys. A: Math. Gen. 21, 2375 (1988). 
[4] J.L. Cardy and I. Peschel, Nucl. Phys. B 300, 377-392 (1988). 
[5] J.L. Jacobsen, J. Phys. A: Math. Theor. 43, 315002 (2010) : ErXiv : 1005 . 36091 
[6] I.G. Enting, J. Phys. A: Math. Gen. 11, 563-568 (1978). 

[7] J.L. Jacobsen and H. Saleur, Nucl. Phys. B 778, 137-166 (2008) ; [arXiv : math-ph/0611078 
[8] J. Kondev and C.L. Henley, Phys. Rev. B 52, 6628-6639 (1995). 

[9] P.W. Kasteleyn and CM. Fortuin, J. Phys. Soc. Jpn. Suppl. 26, 11 (1969); CM. Fortuin and 
P.W. Kasteleyn, Physica 57, 536 (1972). 
[10] H. Saleur, Nucl. Phys. B 360, 219 (1991). 

[11] J.L. Jacobsen and H. Saleur, Nucl. Phys. B 743, 207-248 (2006);[ arXiv: con d-niat /0512058[ 
[12] R.J. Baxter, J. Phys. C 6, L445-L448 (1973). 

[13] F.C Alcaraz, M.N. Barber and M.T. Batchelor, Ann. Phys. 182, 308-343 (1988). 
[14] R.J. Baxter, Proc. R. Soc. Lond. A 383, 43-54 (1982). 

[15] Y. Ikhlef, J.L. Jacobsen and H. Saleur, Nucl. Phys. B 789, 483-524 (2008); 

|arXiv: cond-mat/0612037. 
[16] Y. Ikhlef, J.L. Jacobsen and H. Saleur. larXiv : 11 09 . 1 iW. 
[17] R J Baxter, S B Kelland and F Y Wu, J. Phys. A: Math. Gen. 9, 397 (1976). 
[18] R.J. Baxter, Exactly solved models in statistical mechanics (Academic Press, London, 1982). 
[19] R.J. Baxter, H.N.V. Temperley and S.E. Ashley, Proc. R. Soc. Lond. A 358, 535-559 (1978). 
[20] R.J. Baxter, J. Phys. A: Math. Gen. 19, 2821 (1986). 
[21] R.J. Baxter, J. Phys. A: Math. Gen. 20, 5241 (1987). 

[22] J.L. Jacobsen, J. Salas and A.D. Sokal, J. Stat. Phys. 112, 921 (2003); 

.arXiv : cond-mat/02 04587, 
[23] D. Dei Cont and B. Nienhuis, J. Phys. A: Math. Theor. 37, 3085 (2004); 

[arXiv:cond-mat/0311244l 
[24] J.L. Jacobsen and P. Zinn- Justin, J. Phys. A: Math. Theor. 37, 7213-7225 (2004); 

arXiv : math-ph/0402008, 
[25] D. Dei Cont and B. Nienhuis, [arXiv : cond-mat/0412018| 



Corner free energies and boundary effects 



58 



J. Kondev, Phys. Rev. Lett. 78, 4320-4323 (1997) ; [arXiv : cond-mat/97031 13| 

J.L. Jacobsen and J. Kondev, Nucl. Phys. B 532, 635-688 ( 1998); [arXiv : cond-mat/9804048[ 

J.L. Jacobsen, Conformal field theory applied to loop models. In A.J. Guttmann (ed). Polygons, 

polyominoes and poly cubes, Lect. Notes Phys. 775, 347-424 (2009). 
R.J. Baxter, M.F. Sykes and M.G. Watts, J. Phys. A: Math. Gen. 8, 245-251 (1975). 
I.G. Enting, J. Phys.: Conf. Ser. 42, 83-94 (2006). 

J.L. Jacobsen, J. Salas and A.D. Sokal, J. Stat. Phys. 119, 1153-1281 (2005); 

arXiv : c ond-m at/0401026, 
S. Caracciolo, J.L. Jacobsen, H. Saleur, A.D. Sokal and A. Sportiello, Phys. Rev. Lett. 93, 

080601 (2004); arXiv : cond-mat/0403271, 
T.M. Apostol, Modular functions and Dirichlet series in number theory, 2nd edition (Springer- 

Verlag, New York, 1997), p. 47. 
E. Buffenoir and S. Wallon, J. Phys. A: Math. Gen. 26, 3045-3062 (1996). 
N. Elkies, G. Kuperberg, M. Larsen and J. Propp, J. Alg. Combin. 1, 111-132 (1992); 

arXiv :ma th . CO/920130 5, 
W. Jockusch, J. Propp and P. Shor, l arXiv : math . CO/980 1068[ 

R. Bondesan, J. Dubail, J.L. Jacobsen and H. Saleur, Conformal boundary state for fully open 
boundary conditions, preprint (2011). 



